Witnessing entanglement beyond decoherence 
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We address the witnessing of quantum correlations beyond the limits imposed by an ensemble 
statistical average. By relying upon the continuous observation of a single quantum open system un- 
der the action of classical or quantum noise, we show how the statistical distribution of outcomes of 
physically measurable quantities provides enough information to detect entanglement in situations 
where a simple averaging procedure yields to a vanishing value. In this way, we are able to specify 
the meaning of decoherence when single read-out measurements are available, to demonstrate that 
these distributions are capable to predict /detect entanglement well beyond the ensemble average 
provided by a dephasing-sensitive mixed state. Few-spin and fluorescent systems are used to illus- 
trate entanglement survival when a mixed state average predicts none or under the conditions where 
the sudden-death of entanglement is expected. 

PACS numbers: 03.65.Ud, 03.67.Mn, 75.30.Cr, 42.50.Lc, 42.50.Ar. 



Introduction.- The quantum mechanics most high- 
lighted feature is non-locality, the cornerstone of entan- 
glement, called by Schrodinger "the characteristic trait of 
quantum mechanics" [1]. Locality emerges when a sys- 
tem interacts with an environment (open system), the 
former usually described by a statistical mean, i.e. the 
system's reduced mixed state p, supported by disregard- 
ing beyond the average, the effect from numerous degrees 
of freedom of the latter. As an average, information re- 
garding the spectrum of possible states (decompositions 
or unravelings of p) is degraded, a fact that makes mixed 
state entanglement quantification more complicated and 
less well understood than that of pure states [2]. Re- 
cently, a better understanding of quantum correlations 
in a mixed state has become feasible by addressing the 
entanglement directly linked to physically realizable de- 
compositions Q. Additionally, the increasing inter- 
est to identify measurable observables able to experi- 
mentally witness entanglement (entanglement witnesses, 
EWs) [El, m, jointly with an increased resolution to mon- 
itor the dynamics of single quantum open systems 0, El , 
make it possible to detect quantum non- locality features, 
beyond the smeared out average environment fluctua- 
tions, in optical, atomic/molecular and condensed matter 
physics. 

In this Letter we address the crucial issue of detect- 
ing quantum correlations through the statistical distri- 
bution of EW outcomes captured from continuous envi- 
ronment observation of open quantum systems. Firstly, 
the possibility to witness quantum correlations with de- 
phasing arising from classical noise is studied beyond en- 
semble averages at the level of realizations from a sin- 
gle quantum system, to show that in general the de- 
phasing mechanisms do not always suppress quantum 
correlations, but even more surprisingly, may escalate 
them up. Secondly, an extension to fully quantum en- 
tangled system-environment situations is considered to 
show how witnessing of quantum correlations, beyond 



the mixed-state entanglement-sudden-death (ESD) 
framework, strengthens the relevance of a proper descrip- 
tion of open quantum systems where single read-out mea- 
surements are becoming commonplace [lol. [ll|. 

The mixed state p = ^j^Pk\'^k){'^k\^ is appropriate 
whenever no knowledge on which state \ipk) the sys- 
tem is really in, with infinite decompositions or sets 
{p/e, iV^/c)}, each yielding to a different entanglement 
mean value. A commonly used measure of entanglement 
for bipartite systems in mixed states, the concurrence 
C = Max[0, Ai — XI Ai], where Ai > A^ are square root 
eigenvalues of pp, p = (Jy®(Jyp*(Jy^ (Jy (cr are PauH spin 
matrices throughout this paper), addresses the asymp- 
totic ratio m/n of minimal resource of m ebits required 
to form an ensemble of n ^ oo pairs described by p, and 
singles out the minimally entangled decomposition of p, 
C = min(C) 0,13. For instance, an ensemble of pairs 
of qubits (QBs) with basis states | t), | |) prepared in 
(dephased) fully entangled states (| t|) + e^^| |t))/V2 
with (j) uniformly distributed over (0,27r], has an aver- 
age entanglement (C) = 1 that sharply contrasts with 
the result C = from the purely classical mixed-state 
P = 1(1 tl)(tl I + I lt)Ut I), which arises from averaging 
this ensemble of dephased states. Moreover, any ensem- 
ble averaged physical quantity (therefore, any EW) is, as 
a matter of fact, invariant to any decomposition. Conse- 
quently, it will not unveil the differences in entanglement 
with, say, a half | t|)-half | |t) ensemble ((C) = 0) that 
yields to the same fully mixed state. In order to illustrate 
the relevance of properly relying in EWs statistical dis- 
tributions with realistic physical unravelings from which 
ensemble averages are made upon, spin as well as fluo- 
rescent systems will be addressed below. 

Classical noise.- Consider a system of N interacting 
1/2-spins described by a Hamiltonian {h = 1) such as 
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with distance dependent couplings, specifically dipole- 
like interactions g{rij) = ^o/(^-^x^)^ subject to classi- 
cal noise coming from the random dynamics of the spins 
relative positions rij = \fi — fj\. First, in the absence 
of noise or static case, for = 2 (ri^2 = L) and ini- 
tial state I tii2) (subindex labels each spin), in reso- 
nance UJ2 — oji = A = 0, the concurrence is found to 
oscillate with a frequency go with a full cycle average 
C'mean = 2/7r. On the other hand, far from resonance 
i.e. A ^0, the concurrence oscillates with frequency 
Q = Y^A^ + ^ A reaching a fairly small maximum 
value of Cmax = 2^o^/^^ ^ 2^o/A- Classical noise 
is now included allowing a third particle (3) to elasti- 
cally bounce between the static particles 1 and 2, such 
that ri < (rs) < r2, while it diffuses following a Brown- 
ian motion in an inhomogeneous external magnetic field, 
i.e. co'3(r3) ~ {rs)A/L (Fi^a)). If the translational 
wave packet extension of the mobile particle can be dis- 
regarded, its diffusion is described by a Langevin equa- 
tion {{rs) = rs): r's = — ^r's + ^/wF{t), where F{t) is an 
unit variance white Gaussian noise, r and w are respec- 
tively, relaxation time and diffusion constant, related to 
the diffusivity by D = wr^ /2 [12]. The noisy reservoir, 
here represented by the Brownian particle, leads to de- 
phasing among pure state realizations. The concurrence 
between static spins 1 and 2 is calculated relying upon 
exact numerical simulations for Brownian motion [l2| to 
find the position of particle 3 on a coarse-grained time- 
scale inherent to nanometric particles (r <C This 
last fact implies a single translational time-scale / D 
[16] able to wrap several visits of particle 3 to the static 
ones, useful to enhance entanglement between particles 1 
and 2 when they are out of resonance (Fig(T]^b)). Figure 
[U^c) shows that the noisy motion of particle 3 improves 
both resonance and interaction strength to closely reach 
the noiseless resonant situation and remarkably to en- 
hance over one order of magnitude in average^ the max- 
imum fixed particles original off-resonant entanglement. 
The magnetic susceptibility x is an EW, whenever the 
inequality x/ P < [6] {P is the inverse of temperature 
times Boltzmann constant, 5 = ^ for QBs) holds. A sat- 
uration value of x//^ = | is found for any classical mix- 
ture of states = a\ ^lUU) + 1\ ii^2U) + ^^1 iii2'\?)' 
The closeness of the ensemble average xl P to |, is em- 
phasized in Figdjd), to reflect the fact that the simple 
statistical average (from nearly a classical mixture) is, 
in practice, unable to distinguish the enhancement of 
quantum correlations when noise is present. The maxi- 
mum xl P = \ happens when a = 7 = ^ , and bounds 
a domain down to xl P = | shared with more general 
three QBs classical mixtures. However, the distribution 
of susceptibility from realizations shows (Figdjd)) that 
the main effect of noise is to enlarge the tail of the distri- 
bution on the generally classical forbidden region down 
io xl P = \ (arising from all permutations {a, 7, (5} such 
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FIG. 1: (Golor online) (a) The classical noise physical set up: 
two static particles interacting via dipolar coupling with a 
diffusing third particle (see text), (b) Two typical trajecto- 
ries rzif) (dots) with the corresponding realizations for con- 
currence C (continuous) between spins 1+2. The enhance- 
ment of quantum correlations: the effect of diffusivity D is 
shown in (c), {C)/Cmd.yi reaching closely the average concur- 
rence value for resonant interaction Cmean (arrow); in (d) 
the susceptibility distribution is shown where the value 
xl P — 1.494 zb 0.002 predicted by the average state p for 
D/L^ ^ /\ is highlighted (Qt = 250). In all plots A = lOO^o, 
I '0(0)) = I tii2i3)- In (c) and (d) 5 x 10^ realizations have 
been used. 



quantum correlations is possible from the EW ensem- 
ble statistical average (susceptibility) due to the mixed 
state dephasing-entanglement trade off, is relieved when 
access to the EW statistical distribution is attained, by 
accounting on the additional information available from 
the detailed observation of the environment-induced fluc- 
tuations. 

Quantum noise.- The witnessing of entanglement to 
overcome the inevitable dephasing rooted in ensemble 
averages, can be extended to a fully entangled system- 
environment state 1^), which traced out (hence averaged 
out) over the bath, provides the reduced density operator 
of a mixed state, p = Tr The Schrodinger- 

HJW theorem pj|], states that (infinite) sets of reser- 

1/2 

voir measure-operators {Fj^' } exist without disturbing 
the average state of the system, and whose outcomes 
make visible a specific decomposition of the mixed state 



that a 



7 



hS). The fact that no conclusion on 
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Gonsequently, this last result allows to project onto the 
pure state {i/jk) in which the system is really in, based 
upon the record of environment dynamics. The measure- 
operators produce a change in the system's state when 
the system and environment are entangled, as resolved 
by the so-called Kraus operators ^k{t)^ which fulfill 
Pit) = Ek^iit)p{0)nk{t) and Ek^iitWt) = I. The 
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Kraus operators allow to describe exclusive evolutions 
\ilj{t + At))k = r^/c(At)|V^(t))- preserving purity of states 
from unit efficiency environment measures [1^- condi- 
tioned on a particular outcome k arising with probabil- 
ity P/c = {^k{Aty^k{At))'iij{t)' This evolution (formally 
At is greater than the bath correlation time) proceeds 
with no memory under the Born-Markov approximation, 
according to a Lindblad type master equation 



dtp 



i[H,,p\^Y.^,pL\-\{L\U,p}, 



(2) 



where the operators Li are associated to the irreversible 
effect of the environment on the system. Physically, 
spontaneous emission arising from the interaction of the 
quantized electromagnetic vacuum with an excited atom 
or QB fulfills the Born-Markov approximation. The en- 
semble averaged relaxation from the excited state | t) to 
the ground state | |) of a QB occurring at a rate F is 
described by Eq.(|2]), with L = \/Ta~ . The evolution of a 
single fluorescent QB subject to observation through di- 
rect photo-detection, involves two projectors conditioned 
on no photo-detection Qq = e~i^ t) (t | + | |, 
and photo-detection Cli = Vl — e~^^\ |, imply- 

ing an asymptotic and instantaneous projection onto the 
ground state, respectively [8j. The projectors {flk} can 
be used to construct the Kraus operators for relaxation 
of a pair of independent or interacting QBs (via Hs in 
EqJT]) [1, [HI. In either case, an initial state |?/^o) = tt 
)-\-^\ ii) will present no coherent dynamics. The descrip- 
tion of relaxation subject to photon-detection is provided 
by the Kraus operators Qk=i-\-2j = ^ for indepen- 
dent QBs, while the projection after the first photon- 
detection into unexcited eigenstates \^|^±) = ^(| 
) ± I it)) 7 ill the interacting QBs case, is accounted 
by the replacement Qi^2 ^± = i ^2)- This 

latter decomposition assumes that the quantized vac- 
uum mediated incoherent energy transfer is zero, which 
is exact whenever the induced dipole moments are or- 
thogonal. Interestingly, the no-surveillance mixed state 

Pit) = Er^ritMo){M^Ut) = Es^sitMo){M^lit), 

r{s) = {0, 1(— ), 2(+), 3}, highlights the ambiguous de- 
scription of different quantum states which becomes in- 
variant in an ensambled p level. 

Beyond entanglement sudden death. - The ESD emerges 
when a pair of non-interacting QBs subject to indepen- 
dent dissipative environments, are prepared in a super- 
position |?/^o) that evolves into a mixture with concur- 
rence C{t) = Max[0, e~^^a*(a + e^^{S — a))], vanish- 
ing in a finite time rtssD = ~i^(l ~ whenever 
|<^| ^ 1^1 [9], as shown in Figl2fa). As expected from 
the outlined invariance, the ESD equally occurs with 
QBs interacting via Eq.([T]) (Figl2fb)). However, the use 
of Kraus operators when direct photo-detection is em- 
ployed leads to the abrupt replacement of ESD by an 
asymptotic decay [3] (C)ind = 2|a5*|e~^^ while adding 
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FIG. 2: (Color online) Entanglement of independent QBs (a) 
and interacting QBs (b), initially prepared in the state iV^o), as 
quantified by the analytical expression (C) based on Kraus de- 
composition (thin, black), mixed state concurrence C (thick, 
blue), and intermediate surveillance mean concurrence (Cn) 
from numerical simulations with n = {2, 5, 50} (circles, boxes 
and diamonds, respectively). Inset in (b) shows (C)int ^ind 
C (equal to (C)^-^^) for a = 1. Witnessing of entanglement 
through S"^ is possible as shown from the scatter plot (c), 
becoming an observable able to resolve the distinct entangle- 
ment features from the outcome of on different channels 
for independent QBs (d) and interacting QBs (e) cases, at 
time t — 2tESD. a — such that tmax — tESD- 



a non-monotonic term in the interacting case (C)int = 
2|ape~^^(l — e~^^) + (C)ind to manifest the requirement 
of visiting the maximally entangled superpositions \^|^±) 
that contrasts when a = 1 to the null result of C or 
(C)ind times (Inset Figl2fb)). The usual C from a 

reduced density matrix approach is ill suited to describe 
this case by observing that (C)int presents a maximum 
at Ftmax = ln(2) — ln(l + ■^) only when ESD occurs, i.e. 
if > l^^l- Nevertheless, the comparison among C and 
(C) , reflects the rather abrupt distinct situations of hav- 
ing none or all information possible from the system, due 
to measurement over the environment. 

The effect of having different degrees of available in- 
formation can be addressed if the averaging procedure 
is divided into subensembles of n realizations, pn = 
Xir=i nl^*)(^*l' regarding the mean density operator in 
each subensemble as n ^ 00, while a gradual reduction of 
the subensemble size, until the best resolution is achieved 
when all subensembles have a single realization (n = 1), 
accounts on gaining progressively information concern- 
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ing the times at which individual photons are detected 
[iqI , [i6| . Averaging the concurrence over ah subensembles 
yields to (C^). Figures [2fa)-(b) show how information 
loss due to the averaging of dephased trajectories, swiftly 
replaces the asymptotic value to a finite time concur- 
rence vanishing. However, the robustness against infor- 
mation loss is evident in the independent QBs case, which 
due to averages of factorable states absent of a fragile 
relative phase, present at (C2) no noticeable difference 
with (C) (FiglSfa)), to be compared with the decrease of 
(C2) against (C) in the interacting QBs (Figl2jb)) that 
highlight the delicacy of entanglement quantification due 
to averaging the dephased fully entangled states \i^±). 
These results imply, ahead of conjectures for entangle- 
ment beyond the ESD limit and its enhancement by in- 
coherent relaxation, that quantum correlations in p are 
strongly vulnerable to ensemble averaging, and therefore 
the mixed state can give rise to an unphysical estima- 
tion of the entanglement dynamics in a single quantum 
system subject to noise. 

In order to experimentally verify differences among 
these physical systems, we make use of the photon time 
traces [10| by noting that ESD occurs irrespective of 
spectral overlap with independent QBs while transitions 
to \iIj±) differ by an energy of g{r 1^2)/'^ for interact- 
ing QBs. In this proposal the data from outcomes 
of measuring total angular momentum squared 5^, are 
binned according to the current photon time trace: chan- 
nel (AO) if no photon has been emitted, channel 1 
and 2 (Al and A2), are S'^ measures after the first 
spectrally resolved photon arrival but before channel 3 
(A3), following a second photo-detection. The observable 
S'^ may witness entanglement since its expected value 
(S*^)^ = 2(|ap + |Jp) + + 7p for general pure states 
= a\ tt) +7I ti) + it) + ^^| ii) ranges between two 
and zero, while in fully mixed diagonal states p = a| tt 

)(ni + b\ n)(n 1 + cut)at 1 + (i - a - 6 - c)i u){u 1 

the inequality 1 < {S'^) < 2 holds. Figure [2fc) shows 
that (5^) < 1 only for entangled states, acting in this 
way as an EW, from its calculation on density matrices 
constructed through random normalized complex num- 
bers {a, 7, i^, (5}, with non-diagonal elements further mul- 
tiplied by another complex uniformly distributed number 
within the unit circle such that {i\p\j) = pij = p^^ . 
Hence, the possibility to witness entanglement relies on 
observing {S'^)Ak < 1 in a given subensemble arising due 
to spectral resolution of detected photons. 

Even though the mixed state average (not shown here) 
presents an asymptotical decay from (6'^) = 2 to {S'^) = 1 
unable to verify quantum correlations, its single measure- 
ment read-out will yield to the discrete eigenvalues {0, 2}, 
the former able to witness entanglement. The binned 
discrete outcomes of S'^ for independent QBs (FiglSfd)), 
shows that no conclusion about entanglement can be 
made when S'^ is measured in Al or A2, since in any 
of these both channels {S'^) = 1, regarding the average 



from an equal amount of quantized outputs and 2, ex- 
pected from measurement of total angular momentum in 
a state | t|) or | On the other hand, Figl2fe) shows 
that this measurement is able to witness entanglement in 
the interacting QBs case, since all population in the 
state provides {S'^)a2 = 0. Given that these histograms 
concern a time far beyond the ESD limit, quantification 
of entanglement through C should be disregarded as a 
loyal predictor of entanglement in general decomposi- 
tions. Therefore, this case provides an example where 
witnessing quantum correlations is possible beyond the 
ensemble average when continuous environment observa- 
tion provides the means to classify the outcomes from 
measuring a single quantum system. 

Conclusions.- We have shown how the full statistical 
distribution of physically observable outcomes in single 
open quantum systems, can be used to predict /detect 
entanglement even when the simple ensemble averaged 
EWs, that rely on the dephasing-sensitive mixed state, 
yields to no entanglement at all. The usefulness of single 
read-out EW measurement record from individual quan- 
tum systems, paves the road to experimentally witness 
entanglement beyond smeared out averages. We have 
also shown the need to shift attention from ensemble av- 
erages to physically meaningful decompositions by ad- 
dressing explicitly the worked out examples of a two-spin 
system coupled to a classically diffusing third particle, on 
one hand, and a fiuorescent system, on the other hand. 
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